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Recent work on increasing the effective number of neutrino species (A'^off) in the early universe 
has focussed on introducing extra relativistic species ('dark radiation'). We draw attention to 
another possibility: a new particle of mass < 10 MeV that remains in thermal equilibrium with 
neutrinos until it becomes non-relativistic increases the neutrino temperature relative to the photons. 
We demonstrate that this leads to a value of A'^ofi that is greater than three and that A^cff at 
CMB formation is larger than at BBN. We investigate the constraints on such particles from the 
primordial abundance of helium and deuterium created during BBN and from the CMB power 
spectrum measured by ACT and SPT and find that they are presently relatively unconstrained. We 
forecast the sensitivity of the Planck satellite to this scenario: in addition to dramatically improving 
constraints on the particle mass, in some regions of parameter space it can discriminate between 
the new particle being a real or complex scalar. 



I. INTRODUCTION 

Observations of the abundance of light nuclei pro- 
duced during big-bang nucleosynthesis (BBN) and pre- 
cision measurements of temperature anisotropics in the 
cosmic microwave background (CMB) provide unique 
windows into the particle content of the early uni- 
verse. Both epochs are sensitive to the expansion rate 
of the early universe, which depends on the number of 
relativistic species in thermal equilibrium. The num- 
ber of relativistic species is usually parameterised in 
terms of N^g, the number of effective neutrino species 
with a late-time neutrino-to-photon temperature ratio 
Tv/Try = (4/11)^/^. In the standard cosmological model 
with three neutrino species, this temperature ratio is 
slightly higher than (4/11)^/'^ due to partial reheating 
of the neutrinos when electrons and positrons annihi- 
late, leading to iVcff = 3.046 [T. A measurement of N^g 
greater than this would be a clear sign of new physics. 

BBN imposes limits on iVcft at a photon tempera- 
ture around 1 — 0.1 MeV, since increasing N^g leads 
to a higher abundance of primordial helium and deu- 
terium [2]. While the determination of these light nu- 
clei from astrophysical observations is not yet completely 
settled, the latest data show a preference for a value 
of iVofi greater than three. For instance, Ref. [3] finds 
^cff = 3.7^0 7 (2a) from their inferred abundance of pri- 
mordial "^He. An independent determination of N^s can 
be made at photon decoupling (when the photon tem- 
perature is around 1 eV) from the Silk damping [4] 
tail of the CMB power spectrum. The Atacama Cos- 
mology Telescope (ACT) [5] and the South Pole Tele- 
scope (SPT) [6] have measured this damping tail and. 
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in combination with data from the Wilkinson Microwave 
Anisotropy Probe (WMAP) [7], baryon acoustic oscilla- 
tions (BAG) [8] and the Hubble constant (Fq) PHH], also 
find a preference for a value of iVoff greater than three, 
obtaining iVeff = 4.6 ± 0.8 (Ict) and N^g = 3.9 ± 0.4 (Ict) 
respectively. Furthermore, various independent analyses 
of recent cosmological data find evidence for Nqq > 3 at 
95% CL [IIHISI- 

Given these hints for extra energy density in the 
early universe and bearing in mind the expectation that 
the Planck experiment will soon provide a significantly 
improved measurement of N^s at photon decoupling 
(e.g. see |14j). we consider it a pertinent time to con- 
sider models in which iVcff is increased. The conven- 
tional way is to introduce extra 'dark' radiation, which 
leads to the same value of iVoff at BBN and photon 
decoupling (see e.g. [TSHlT] ). While it is usually as- 
sumed that A'cff does not change between BBN and pho- 
ton decoupling, intriguingly, the experimental data are 
consistent with a slightly larger value of A'cff at pho- 
ton decoupling. This increase in Ni^s after BBN can be 
achieved through decays to dark radiation during or after 
BBN (see e.g. [THH15]) or from primordial gravitational 
waves (211 . 

In this paper we explore another scenario that has re- 
ceived less attention in the literature, which predicts a 
value of iVeff greater than three, and moreover, predicts 
that iVeff at photon decoupling is larger than the value 
during BBN. As we will show in Section [llj TVeff scales as 
(Ty/T-y)*^ so increasing the neutrino-to-photon tempera- 
ture ratio leads to an increase in A'cff. A relative increase 
in T^/T^ can be achieved by reducing the photon temper- 
ature relative to the neutrino temperature or by increas- 
ing the neutrino temperature relative to the photon tem- 
perature. A decrease in can be obtained through the 
production of hidden photons [251 ES] while a new light 
mediator may lead to a change in either or T^, ^7] . In 
this paper we follow the example of the standard cosmo- 
logical model to increase : in the standard cosmological 
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model, the photons are hotter than the neutrinos because 
the transfer of entropy from the electrons and positrons 
to the photons (when they become non-relativistic) hap- 
pens after the neutrinos decouple from the electromag- 
netic plasma at To « 2.3 MeV [55]. In an analogous 
way, we show that an additional 'generic' particle x that 
remains in thermal equilibrium solely with the neutrinos 
after they decouple from the electromagnetic plasma and 
until it is non-relativistic, reheats the neutrinos relative 
to the electromagnetic plasma and as a result, leads to 
a higher value for N^s- Requiring that the neutrino re- 
heating, which happens when x becomes non-relativistic 
and transfers its entropy to the neutrinos, happens af- 
ter neutrino decoupling implies that x must have a mass 
ruy < few • Td ~ 10 MeV. 



The goal of this paper is to investigate the current and 
future constraints on x from the inferred values of N^s 
from BBN and the CMB . The only condition that we 
impose on the 'generic' particle x is that the neutrino- 
X interaction rate is sufficiently high that they remain in 
thermal equilibrium until x is non-relativistic. Therefore, 
the constraints we derive may apply to, for instance, light 
dark matter particles that obtain their abundance from 
thermal freeze out or light mediators that have large cou- 
plings to neutrinos. The parameter that determines N^s 
is as this dictates the additional energy density, so it 
is this parameter that we will constrain. During BBN, 
the photon temperature is similar to so x makes a di- 
rect contribution to A^off • In addition, there may also be 
an indirect contribution from the increase in Ti,/T^. At 
photon decoupling, we assume x is non-relativistic so its 
direct contribution to 7V(,g is Boltzmann suppressed; the 
increase in N^a arises solely from the increase in TyjT^. 
As we show in Section |llj this difference in the origin of 
the extra energy density leads to a larger value of iVoff at 
photon decoupling. 



This paper is organised as follows. In Section|n]we dis- 
cuss the impact of light particles in thermal equilibrium 
with neutrinos on the energy density and neutrino-to- 
photon temperature ratio (some results are fully derived 
in an Appendix). Using these results, we first find the ef- 
fect of X on the abundance of primordial nuclei produced 
during BBN. This has previously been studied in [29, 30J. 
Here, we independently calculate the primordial abun- 
dance of ^He and D/H as a function of and compare 
with recent experimental measurements. Our calculation 
uses a more recent BBN code than j29l [30] and includes 
the most recent values of the neutron lifetime and baryon 
density. We then calculate the value of A^eff as a func- 
tion of at photon decoupling and compare with the 
experimental values inferred by ACT and SPT. Turning 
to consider experimental results from the near future, in 
Section UTTl we forecast the constraints that will soon be 
placed on x from Planck's measurement of iVeg. Finally, 
we conclude in Section ITVl 



II. CURRENT CONSTRAINTS 

The total energy density of the universe /Or is usually 
parameterised in terms of the energy density of photons 
p-y, and the effective number of neutrinos f^cS with the 
usual neutrino-to-photon temperature ratio jT^ (we 
define this ratio more carefully below) 



PR = Pi 



7 /T" 

1 + - — 
8 



off 



(1) 



We are interested in the case where an additional par- 
ticle X with internal degrees of freedom and mass 
is in thermal equilibrium with the neutrinos. As the par- 
ticles are in thermal equilibrium, they share a common 
temperature T^, which will in general be different from 
T^. The resulting energy density of N^, neutrinos and 
the additional particle x is 
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where for convenience, we have defined y^, = m-^/T^ and 
the function 
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which takes the limits I{y — > oo) = and 
I{y — )■ 0) = 1 (8/7) for fermions (bosons) (as usual, in 
Eq. ([3|, —1 pertains to bosons and -1-1 to fermions). 
Comparing Eq. ([2| with the definition of A'off in Eq. ([I]) , 
we see that 
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In the Appendix, we use the conservation of entropy 
per-comoving- volume for particles in thermal equilibrium 
to show that 
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where Qi^s-.u and g-ks-.-i are the effective number of rela- 
tivistic degrees of freedom in the neutrino and electro- 
magnetic sectors respectively and \rp^ indicates that g-i,s 
should be evaluated at the neutrino decoupling tempera- 
ture Td- In the absence of x, g*s:v remains constant after 
decoupling, hence the usual neutrino-to-photon ratio is 
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More generally, we see that 
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If X becomes non-relativistic after decoupling, 
<?*s:i/ < g*s-M\T-a and as a result, the neutrino-to- 
photon temperature ratio increases above its usual 
value. From Eq. Q, we see that this leads to an increase 
in A^^eff- 

Finally, in the Appendix we show that Eq. ([7| can be 
written as 
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where Hi^lxj^ = m-^/Ti) and we have defined the function 
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which takes the limits F{y — > oo) — and 
F{y — >■ 0) = 1 (8/7) for fermions (bosons). It is this form 
of Tu/T^ that we use in our calculations. In the following 
subsections, we use these formulae to explore the impli- 
cations of X on BBN and the CMB. 



A. Effect on BBN 

The abundance of "^Hc, expressed through the mass 
fraction Yp, has long been recognised as a probe of the 
energy density present during BBN |3T] . In the early uni- 
verse the ratio of the neutron-to-proton number densities 
(n„/np) is kept in thermal equilibrium through the weak 
interactions until the expansion rate becomes comparable 
to the weak interaction rate at T ^ 0.7 MeV. Increasing 
the expansion rate, by for instance introducing additional 
relativistic radiation, leads to a larger value of n„/np at 
freeze out. Since essentially all of the neutrons present at 
T ~ 0.1 MeV are synthesised into ^He, introducing ad- 
ditional relativistic radiation therefore leads to a larger 
abundance of "^He. Similarly, increasing the expansion 
rate means that the reactions depleting the abundance of 
D freeze out earlier, leading to an increased abundance. 
The reactions that deplete D are more sensitive to the 
baryon-to-photon number density 77 and have historically 
been used to infer its value. However, if we use the value 
of 77 inferred by WMAP, then the number ratio D /H can 
be used along with 1^ as a measure of the energy density 
present during BBN, a point emphasised in [T3] . Finally, 
the observed abundance of ^Li remains ~ 5(t away from 
the theoretical prediction, the so-called 'lithium problem' 
(for a recent review, see 32J). The modifications that 
we propose do not significantly alter the ^Li abundance, 
therefore, in this work we directly compare our theoreti- 
cal predictions with the inferred values Yp and D /H from 
astrophysical measurements and make only passing ref- 
erence to ^Li. 

The impact on the primordial abundances of light 
nuclei from a generic particle x that maintains ther- 
mal equilibrium with neutrinos throughout BBN has 



previously been considered in [551 ISO]- Here we in- 
dependently calculate Yp and D/H using a modi- 
fied version of the PArthENoPE BBN code [33]. We 
use the latest measurements of the baryon density 
Vl\jh? = 0.0223, deduced from a joint analysis of 
SPT-|-WMAP7-h7?o+BAO j6|, and the PDG value of the 
neutron hfetime r„ = 880.1 s [Mj . 

In the spirit of |30j . we introduce a small temperature 
dependent parameter 



5{T) 



1 



r°(r) 



(10) 



and perturb the PArthENoPE code in order to take into 
account the effects from the modified neutrino-to-photon 
temperature ratio. We use Eq. (|8| to calculate T°/T^. 
Typically, 5 ^ 0.01 and is always smaller than ^0.1. A 
non-zero value of 5 increases the energy density of the 
neutrinos and enters into the phase-space of the weak in- 
teraction rates that determine n„/np. The PArthENoPE 
interaction rates r„_i.p and Tp^n include finite mass, 
QED radiative and finite temperature corrections. We 
calculate corrections to these rates in the Born approxi- 
mation as an expansion in 5 such that the total rate for 
n — >■ p is given by 
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where eni{T) fits the change in the rates to an accuracy 
better than a few percent. A similar expression holds 
for the p ^ n rate. Finally, we also include the extra 
contribution to the energy density from x, expressed in 
Eq. ([2|. 

The red dot-dashed, solid and dashed lines in Fig. [T] 
show our results for a complex scalar (B2), Majorana 
fermion (F2) and real scalar (Bl) respectively as a func- 
tion of m^. The upper and lower panels show the values 
of Yp and D/H respectively and are in good qualitative 
agreement with those found in [221 [30] , with slight differ- 
ences due to the updated parameter values that we use. 
Although not shown here, our prediction for ^Li is sim- 
ilar to that in [30]. For reference, the black dotted lines 
show the predicted values of Yp and D/H for the indi- 
cated values of Ncs- As we would expect, for > 15 
MeV, we recover the result from standard BBN: this is 
because x is non-relativistic during BBN (so its contribu- 
tion to the energy density is Boltzmann suppressed) and 
it has transferred its entropy to the neutrinos before they 
decouple from the photons, meaning that the standard 
neutrino-to-photon temperature ratio relation is main- 



tained {gi,s-M = g*s:f\Tj, in Eq. Q). For < 0.05 MeV, 
we asymptote to the result expected from a massless par- 
ticle as X remains relativistic throughout all of BBN. The 
values of Yp and D/H for intermediate values of are 
a result of the direct contribution to the energy density 
from X and the contribution from the modified neutrino- 
to-photon temperature ratio. 

Recent inferences of Yp from observations of metal- 
poor H II regions have been slightly higher than re- 
sults from the past decade. For instance, while Ref. [35] 
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FIG. 1: The red dot-dashed, solid and dashed lines show the 
predictions for Yp (upper panel) and D/H (lower panel) for a 
complex scalar (B2), Majorana fermion (F2) and real scalar 
(Bl) respectively. The blue shaded region indicates the lo 
region for Yp from [37] and the la weighted mean of D/H 
from [39] . The black dotted lines show the values of Yp and 
D/H for the indicated values of A'j.ff. 
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FIG. 2: The red dot-dashed, solid and dashed lines show A'efi 
at photon decoupling ( A'eff Iqi^jb) as a function of m-^ for a 
complex scalar (B2), Majorana fermion (F2) and real scalar 
respectively. The shaded green and blue hatched regions show 
the Icr range of A^oft determined from ACT and SPT (in com- 
bination with data from WMAP7, BAO and Hq) respectively. 



ture, a recent preprint provides a very precise measure- 
ment of D/H with an error ~ 4 times smaller than the 
previous best |42j . With this measurement, they find 
that the weighted mean from ten independent samples is 
D/H X 10^ = 2.63 ± 0.12 (where the shift in the weighted 
mean and reduction of the error is largely due to the new 
measurement in [42j ) . Comparing with the lower panel of 
Fig. [ij we see that this result disfavours < 4 MeV at 
2(7 for a Majorana fermion and complex scalar while still 
allowing all values of for a real scalar. With further 
observations from metal-poor damped Lyman alpha sys- 
tems of the type considered in (32 , the large dispersion 
in the measurements of D /H would decrease and provide 
a strong, reliable constraint on this scenario. 



found Yp = 0.249 ± 0.009, more recently, Refs. [3], 
[5S], [37] and [35] found Yp = 0.2565 ± 0.0010(stat.) ± 
0.0050(syst.), Yp = 0.2561 ± 0.0108, Yp = 0.2573^g;S5!]^| 
and Yp = 0.2534 ± 0.0083 respectively, in good agree- 
ment with each other. As a representative of the more 
recent measurements of Yp, the blue shaded region in the 
upper panel of Fig.[I]shows the Icr region for Yp from [57] . 

The experimental value for the number ratio D/H 
shows a large spread with the dispersion in the measure- 
ments well in excess of the quoted errors. For example, 
measurements range from D/H x 10^ = 1.17lg'34 in [ID] 



to D/H X 10^ 
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Here we take as a 



conservative value D/H x 10^ = 2.78lo:2i, the weighted 
mean calculated from eight independent samples in [55] . 
The blue shaded region in the lower panel of Fig.[l]shows 
this la region. It is clear from Fig. [l] that the error on 
the measurements of Yp and D/H are sufficiently large 
that all values of are currently consistent with the 
data at just over la. Although the measured errors of 
Yp are unlikely to improve considerably in the near fu- 



B. Effect on the CMB 

An independent determination of can be made 
from the CMB, which is at a much lower photon tem- 
perature than found during BBN. For the values of 
that we consider, x's direct contribution to the energy 
density is Boltzmann suppressed at photon decoupling 
(T.7-dcc 1 eV): that is, I{m^/T^.Aec) ~ in Eq. ([2|. 
However, assuming that x remains in thermal equilib- 
rium with neutrinos until it is non-relativistic, we can 
calculate its indirect contribution from the increase in 
the neutrino-to-photon temperature ratio from Eq. ([8| 
and the resulting value of NcS from Eq. At photon 
decoupling we have F{m^/T^.iicc) — so A'off is given by 
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where we have set = 3.046 [1. In Fig. |2] we show 
iVoff as a function of for a complex scalar (B2) , Majo- 
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rana fermion (F2) and real scalar (Bl) in red dot-dashed, 
solid and dashed lines respectively. For ^15 MeV 
we recover the usual value for Ncs because the reheating 
effect from x takes place before Td , so the neutrinos and 
photons are both reheated equally. For m-^ < 2 MeV, 
A^off remains constant because the neutrinos are maxi- 
mally reheated for all of these values of m^. This is 
because x becomes non-relativistic after Td and is fully 
non-relativistic by photon decoupling, therefore, all en- 
tropy from X is transferred to the neutrinos. Naively, we 
might have expected A^cff to asymptote to four when x 
is a light Majorana fermion. That this does not happen 
can be traced to the fact that the term in square brackets 



in Eq. (12) is raised to a power 4/3, rather than 1. This 
in turn follows from entropy conservation; the factor of 3 
follows from the fact that S ^ gi,sT^ is conserved, rather 
than g^,sT^. 

Experimentally, -^cfflcMB determined from the 
CMB power spectrum at high multipole £ ('the Silk 
damping tail'). Increasing A^off increases the energy 
density, which in turn increases the expansion rate. 
As discussed in (43) . an increased expansion rate leads 
to increased Silk damping which reduces the power in 
damping tail of the CMB power spectra. Currently, 
ACT [5] and SPT [6 provide the most precise mea- 
surements of the CMB power spectra at high t (up to 
I ^ 3000) and are thus able to set the best limits on 
-^offlcMB from the CMB power spectrum. In Fig. [2j 
we show the current la constraints on Nefi\cMB, froni 
ACT (shaded green region) and SPT (blue hatched re- 
gion) when their respective data sets are combined with 
data from ^MMKPI+BAO+Hq} While we see that the 
ACT and SPT data show a preference for N^s > 3 and 
therefore, < 10 MeV, all values of are currently 
consistent with the data at 2cr. As a result, no strong 
conclusions on the preferred values of can currently 
be drawn from the data. 
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FIG. 3: The difference in A'eff at photon decoupling and BBN 
(inferred from the value of Yp). The red dot-dashed, solid and 
dashed lines show the prediction for a complex scalar (B2), 
Majorana fermion (F2) and real scalar (Bl) respectively. 



mains in thermal equilibrium with neutrinos until non- 
relativistic, we see that -/VofrlcMB > A'off [l^JIggj^. The 
equality holds for 15 MeV because A'cff remains at 

its usual value. For < 0.2 MeV, the difference is es- 
sentially constant because A^eff takes the maximum value 
at BBN and photon decoupling. The difference in TVoff is 
greatest for intermediate values of . 

Finally, for completeness we mention a way to achieve 
A^cff IcMB ^ ^off [ypllBBN- If X decouples from the neu- 
trinos while still relativistic but after BBN (implying 



< 0.1 MeV) we would have A'og [I^IIbbj^ — 4 and no 
reheating of the neutrinos would occur. Furthermore, if 
> 1 eV, X would not contribute to the energy density 
at photon decoupling since it would be non-relativistic. 
Therefore, -/VosIcmb — 3.046 as usual. We will not ex- 
plore this possibility any further. 



C. Comparison 

Comparing Figs, [l] and [2j we see that A^cfflBBN ^^"^ 
CMB differ for the same value of m^. In Fig. p^we 
explicitly show the difference in N^s between the BBN 
and the CMB epochs as a function of for a complex 
scalar (B2), Majorana fermion (F2) and real scalar (Bl) 
(in red dot-dashed, solid and dashed lines respectively). 
Here, we follow the usual convention and use the inferred 
value of A^offlgBfj from Yp (-/Vcff [^]Ibbn) but note that 
a similar figure could be drawn in which A'offlBBN i^ i'^" 
ferred from the values of D/H. When a particle x ^e- 



Strictly speaking, these constraints apply when Yp is fixed to the 
standard BBN value with N^s massless fermions. However, from 
Table [I] we see that the constraints on Ncs\cMS (from SPT) do 
not significantly change when Yp is left as a free parameter. 



III. PLANCK FORECAST 

In the previous section, we demonstrated in Figs, [l] 
and [2] that all values of are consistent (within 2a) 
with the current measurements of A'off from BBN and 
the CMB for a particle Xi which remains in thermal equi- 
librium with the neutrinos after they decouple from the 
electromagnetic plasma and until it is non-relativistic. 
However, we would ideally like to place strong constraints 
on or find evidence for such a scenario. As Planck will 
soon provide significantly improved measurements of the 
CMB temperature anisotropy, we now turn to forecast 
its sensitivity to our scenario. 

In Section [il B| we stated that A^cff is determined from 
measurements of the damping tail of the CMB power 
spectrum. Due to the accuracy of Planck, we must be 
careful to take into account other parameters that also 
lead to additional damping of power on small scales as 
they will be degenerate with N^s [44]. For instance, a 
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negative running of the spectral index Ug with scale k 
is one such parameter. However, the value required to 
provide sufficient damping, dris/dlnk ~ —0.02, is much 
larger than expected from slow-roll inflationary models 
and would require a rethink of inflation and the genera- 
tion of perturbations [551 ES] so we will not consider it 
further. Other examples of such parameters like an alter- 
native model of dark energy, are considered in [47l [48] . 
Here, we concentrate on a conventional, flat, ACDM 
model that we extend to include a parameter that has 
already been central to our discussion, namely, Yp. Since 
*He recombines earlier than hydrogen, increasing the *He 
abundance (at fixed baryon density) leads to fewer free 
electrons during hydrogen recombination. This in turn 
leads to a larger photon diffusion length and as a result, 
less power in the CMB damping tail |49] . Increasing N^s 
can be compensated by decreasing Yp so these parameters 
are degenerate [il]. 

In the first column of Table |l] we list the cosmologi- 
cal parameters that we specify in our forecasts. In ad- 
dition to rib^^, TT-s, Ncff and Yp which have previously 
been defined, we choose values for the cold dark matter 
density ^luuh'^', the reduced Hubble parameter h; the 
angular scale of the sound horizon at last scattering 9s f 
the optical depth to reionization t; the amplitude of pri- 
mordial scalar fluctuations ln[10^''As] (at a pivot scale 
k — 0.002 Mpc^^) and the fraction of dark matter in the 
form of neutrinos = ri^/f^DM- Finally, we also list the 
parameter erg, the amplitude of linear matter fluctuations 
on scales of 8h~^ Mpc at z = 0, which is constrained by 
galaxy clusters. 

In order to find the current best fit values for our 
parameters, we perform a fit with CosmoMC [50] to the 
SPT-^WMAP7+BAO+ifo data. We use v4.1 of the 
WMAP likelihood code and vl.2 of the SPT likelihood 
code publicly available from the LAMBDA website \5T\ . 
Following when calculating the SPT likelihood we 
also sample over three foreground nuisance parameters, 
in addition to the cosmological parameters listed in Ta- 
ble U The results of our fit are shown in the second col- 
umn of Table H 

To generate mock CMB data, we follow [511 [53] and 
use the FUTURCMB [Si] add-on package to CosmoMC. We 
use the CAMB software [5F package to obtain angular 
power spectra for the CMB anisotropics. From this we 
use the TT, TE and EE angular spectra for multiples 
2 < £ < 2500, and the lensing deflection maps for dd 
and Td. Planck's i3-mode measurement will likely be 
noise-dominated so we omit it from our fits. We assume 
a sky coverage of /sky = 0.65 and follow in our ap- 
proach to forecasting by replacing the power spectrum of 
the mock data by the fiducial data set. Using CosmoMC, 



^ In generating the mock data with CAMB, we choose h and 9s 
is derived. In exploring the parameter space with CosmoMC, we 
sample over 6s and h is derived. For completeness we hst the 
values for both parameters. 
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U.O — T 






n I'^ri-i-nmiQ 
u.iou It u.uiiy 


0.111 


0.134 


n m 
u.ui — 


^ n QQ 


n 




0.71 


0.76 


U.4 — ? 


1 n 

i.U 


100 6s 


1.040 ± 0.0030 


1.043 


1.039 


0.5 -> 


10 


■J- 


0.090 ± 0.015 


0.09 


0.09 


0.01 - 


> 0.8 




3.17 ±0.04 


3.16 


3.13 


2.7^ 


4 


ns 


0.989 ±0.013 


0.98 


0.99 


0.5 ^ 


1.5 




< 0.054 


0.008 


0.020 


^ 1 




iVcff 


3.98 ±0.61 


3.046 


4.418 


2 -> 7 






0.2629 ±0.039 


0.247 


0.257 


0.22 - 


> 0.27 




0.801 ±0.054 


0.803 


0.816 







TABLE I: The first column lists the cosmological parameters 
we sample over in our Planck forecast. The second column 
shows the values we obtain for these parameters when fitting 
to the SPT+WMAP7±BAO±Ji'o dataset. For two-tailed dis- 
tributions the 68% CL is given while for one-tailed distribu- 
tions, the upper 95% CL is given. In the third and fourth 
column we list the fiducial values of the various parameters 
chosen to generate the mock CMB data. 'ACDM' corresponds 
to a cosmology without x while in 'ACDM±x', and A^eff 
take the values for a 1 MeV Majorana fermion. The last col- 
umn lists the flat priors we use in CosmoMC. 



we generate eight Markov chains in parallel and monitor 
the convergence with the Gelman- Rubin _R-statistic |56j . 
requiring that i? — 1 < 0.015. 

In the third and fourth columns of Table H] we list the 
fiducial values of the various parameters that we use to 
generate our mock CMB data. These are chosen to be 
consistent at or just outside the Icr error range of the 
parameters in the second column. In the third column 
('ACDM'), Yp and A^eff are set to the values predicted 
from the standard concordance model. In the fourth col- 
umn ('ACDM+x'), Yp and A'off take the values expected 
from a 1 MeV Majorana fermion (as determined from 
Figs, [l] and [2] respectively) . In the fifth column we list 
the flat priors which we impose on the various parame- 
ters when exploring the parameter space with CosmoMC. 
Since ./Veffjggj^ ^ -^osIcmB' longer have the usual 

relation between Yp and N^n so we do not impose a BBN 
prior on Y^, as suggested in [S3]. Rather, we choose a 
conservative flat prior on Yp, determined from the cen- 
tral value of Ref. [37] {Yp = 0.2573) ±4cr. Finally, 
we have chosen parameters to be within 2a of the ob- 
served value of as: Ref. [57] found cr8(f^m/0.25)°-4^ = 
0.813 ± 0.013(stat.) ± 0.024(syst.), where is the mat- 
ter density. We have that CT8(f^m/0.25)°-4'^ = 0.827 and 
0.855 for 'ACDM' and 'ACDM+x' respectively. 

In Fig. [4] we show the results of our forecast in the 
Yp-NcB plane. The blue contours in the upper and lower 
panels of Fig. [4] show the la and 2a regions we find for 
Planck, after sampling over the mock power spectra cre- 
ated with the 'ACDM' and 'ACDAI+x' fiducial values 
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^effl 



FIG. 4: The top and bottom panels correspond to the fidu- 
cial values of 'ACDM' and 'ACDM+x' from Table |l] respec- 
tively. Both panels: In blue are the Planck la and 2a re- 
gions. The red lines show the values of 1^ as a function of 
A'^eft for a real scalar (Bl), complex scalar (B2) and Majorana 
fermion (F2). The black star indicates the values of Yp and 
A^cff for the stated values of in the box. Upper panel: In 
green are the SPT-|-WMAP7-|-BAO-|-//o lo" and 2a regions. 
Lower panel: In green are the la and 2a regions for Yp from 
Aver et. al. ^7]- The black dotted line marked DR indicates 
the relation between Yp and N^b for dark radiation in which 

A^ohIbbN = ^EfflcMB- 



from Table |T] respectively.^ In both panels, the red dot- 
dashed, solid and dashed hnes show the relation between 
Yp and A'eff in the Yp-N^g plane for a particle x when 
it is a complex scalar (B2), real scalar (Bl) or Majorana 
fermion (F2) respectively. Comparing with Figs. [T] and 
[2] it is clear that the left (right) edge of the lines corre- 
spond to the values of Yp and A'eff for large (small) values 
of m^. 



^ We do not show the one-dimensional distributions for the in- 
dividual parameters as we are primarily interested in how well 
Planck can discriminate between different values of . The one 
dimensional distributions we find have a similar variance to those 
found in, for instance [53) . 



In the upper panel of Fig. [4] in green we 
show the Icr and 2cr regions after fitting to 
SPT-HWMAP7+BAO+i7o- Confirming the conclusions 
from Section|TlB] we see that SPT-t-WMAP7+BAO+i7o 
does not place strong constraints on m^, with all values 
lying within the 2cr region. In comparison, Planck's 
confidence regions are much smaller and allow much 
more discrimination between values of m^. For instance, 
if Planck makes a measurement of A'eff close to three, 
as in the upper panel of Fig. |4j we see that it is able 
to exclude at 2a values of below 4.8 MeV for a real 
scalar and 8.0 MeV for a complex scalar and Majorana 
fermion. 

In the lower panel of Fig.|4j the green regions show the 
1(7 and 2a regions from the inferred value of Yp from Aver 
et. al. [37]. The black dotted line marked DR indicates 
the usual relation between Yp and A'eff for dark radia- 
tion in which A'eff Iggp^ = A'eff | ^j^g . The deviation from 
this relation for a real scalar, complex scalar and Majo- 
rana fermion occurs because A'eff Ipfyig > A'eff [Yp] Iggjj. 
If Planck finds evidence for a larger value of A'eff , we see 
that it will be challenging to distinguish between the case 
of dark radiation and a generic particle x with high sta- 
tistical significance unless the value of Yp inferred from 
astrophysical measurements improves significantly. How- 
ever, considering only particles in thermal equilibrium 
with neutrinos, we see that Planck is able to discriminate 
between a real scalar and a complex scalar or Majorana 
fermion at more than 2a and is able to exclude at 2a 
values of heavier than 4.2 MeV for a complex scalar 
and 3.8 MeV for a Majorana fermion. 



IV. CONCLUSIONS 

Motivated by discrepancies in the determination of 
Acff from BBN and the CMB, we have considered the 
impact of a new light particle which remains in ther- 
mal equilibrium with neutrinos until it becomes non- 
relativistic, which we assume occurs before photon de- 
coupling. Such a particle leads to extra energy density 
in the early universe, either directly or through its effect 
on the ratio of the neutrino-to-photon temperature. 

To quantify the impact of such particles on BBN, we 
updated the analysis of [29j [30] through modifications 
of the PArthENoPE code to take into account the two 
effects of light particles on the energy density mentioned 
above. These particles lead to an increased value of A'eff 
and bring the ^He abundance into better agreement with 
recent observations (see the upper panel of Fig. [l]). We 
also showed (in the lower panel of Fig. [T]) that the D 
abundance is compatible with recent measurements. 

We then considered the effect of such particles at the 
time of formation of the CMB. At this epoch, the increase 
in A'eff comes solely from the increase in the neutrino-to- 
photon temperature ratio, as x is non-relativistic. As 
we showed in Fig. [2] this brings the value of A'eff into 
better agreement with the values reported by CMB ex- 
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perimcnts. In general, we found that the value of N^g at 
the time of CMB formation is larger than the value at 
BBN. This is demonstrated explicitly in Fig. [3] Fixing Yp 
to the value inferred from iVoff at photon decoupling, as 
is often done, is not applicable when setting constraints 
on this scenario. 

As the current experimental constraints on our sce- 
nario are still relatively weak, we forecast the sensitivity 
of the Planck satellite to the effect of these particles on 
the CMB. We find that Planck is highly sensitive to the 
effects of such particles. If the value of iVcff derived from 
Planck agrees with that of standard cosmology with three 
neutrinos, we showed in the upper panel of Fig. |4]that it 
will rule out (for example) the existence of a light Majo- 
rana fermion with a mass of less than 8 MeV which cou- 
ples to the neutrinos. The lower panel of Fig. |4] indicates 
that it will be difficult to distinguish between our sce- 
nario and the standard case of dark radiation (in which 
A'eff is the same at BBN and photon decoupling) at much 
more than la when only data from astrophysical mea- 
surements of ''He and Planck is considered. On the other 
hand, considering only particles in thermal equilibrium 
with neutrinos, we demonstrated in the lower panel of 
Fig. |4] that as well as providing significantly improved 
constraints on m^, in some regions of parameter space, 
Planck is even able to discriminate between a real scalar 
and a complex scalar or Majorana fermion. 
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with common temperature T. That is 



S 



3 P + P 
T 



constant 



(A.3) 



Here, a is the scale factor and this equation defines 5*5, 
the effective number of relativistic degrees of freedom. 

At r = Td , the neutrinos and x (the neutrino plasma) 
decouple from the electrons and photons (the electromag- 
netic plasma). From this point onwards, the neutrino 
and electromagnetic plasmas separately conserve comov- 
ing entropy. As a result, the neutrino and photon tem- 
peratures (T^, and respectively) evolve according to 



rp _ 1 -1/3 _1 



(A.4) 
(A.5) 



where ki and k2 are constants and gi,a:i> and g^si^ are the 
effective number of relativistic degrees of freedom in the 
neutrino and electromagnetic plasmas respectively. As- 
suming that the neutrinos and Xj a-nd the photons and 
electrons decouple instantly at = = Td (and ignor- 
ing that muon and tau neutrinos decouple slightly earlier 
than the electron neutrino [SH]), we can determine fci/fc2 
and therefore, we find (for < Td) 



9*s:f 




(A.6) 



where \rp^ indicates that g^^s should be evaluated at Td. 

Using Eqs. ( A.l ) and ( |A.2[ ), the entropy per comoving 
volume for N^, neutrinos and x is 



3 27r^ 3 
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(A.7) 



Appendix: Neutrino-to-photon temperature ratio 

The energy density p and pressure p of a particle of 
mass m at temperature T and g internal degrees of free- 
dom are given by (ignoring a possible chemical potential) 



g_T^ 
27r2 



P - 



'4 I'oo 



e« ± 1 



67r2 Jy ^ e« ± 1 



(A.l) 
(A.2) 



where ^ = E/T, y = m/T and here, as usual, —1 pertains 
to bosons and -1-1 to fermions. 

In order to calculate the neutrino-to-photon temper- 
ature ratio, we use the conservation of entropy per- 
comoving-volume S for particles in thermal equilibrium 



where yi, = m^/Ti, and for convenience, we have defined 
the function 



30 
7^ 



e« ± 1 



(A., 



with limits, F{y ^- 00) = and F{y ^ 0) 1 (8/7) for 
fermions (bosons). Comparing Eq. (A.7) with Eq. ( A.3[ ), 
it is straightforward to read off that 



g*s:u — 2 • 



A similar argument shows that 



7 

g^s:^ = 2 + 4 • - • F(z, 



(A.9) 



(A.IO) 



where = me/Ty 
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